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ABSTRACT 

Using the Keck I telescope, we have obtained 11.7 /im and 18.7 /im images of the circumstellar dust 
emission from AB Aur, a Herbig Ae star. We find that AB Aur is probably resolved at 18.7 //m with 
an angular diameter of 1.2" at a surface brightness of 3.5 Jy arcsec -2 . Most of the dust mass detected 
at millimeter wavelengths does not contribute to the 18.7 /im emission, which is plausibly explained 
if the system possesses a relatively cold, massive disk. We find that models with an optically thick, 
geometrically thin disk, surrounded by an optically thin spherical envelope fit the data somewhat better 
than flared disk models. 

Subject headings: stars: individual (AB Aur) — circumstellar matter — planetary systems: formation 
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1. INTRODUCTION 

Herbig Ae/Be stars are pre-main sequence progenitors 
of intermediate mass A- and B-type stars which possess 
strong infrared excesses, dominated by thermal emission 
from dust grains. It is generally thought that these stars 
possess circumstellar disks (Mannings & Sargent 1997, 
2001) in which planets might form. The grain properties 
(size, temperature and composition) and the macroscopic 
geometry of the mid-infrared emitting dust around these 
stars are not yet well determined. In the simplest mod- 
els, Herbig Ae/Be stars possess optically thick, geomet- 
rically thin disks which produce spectral energy distribu- 
tions (SEDs) with vF v oc z^ 4 / 3 , steeper than the observed 
SEDs (Hillenbrand et al. 1992). Several groups have sug- 
gested a variety of explanations for this discrepancy such 
as flared disks with a surface layer of small grains (Chiang 
& Goldreich 1997) and geometrically thin disk embedded 
in spherically symmetric envelopes (Miroshnichcnko ct al. 
1999). We have carried out a high resolution mid-infrared 
imaging study of AB Aur to learn about the spatial dis- 
tribution of warm dust which may allow us to distinguish 
among models. 

AB Aur is a AO Ve+sh (shell) Herbig Ae star with a 
Hipparcos distance of 144 pc from the Sun. Since AB Aur 
lies near the zero age main sequence, there is some uncer- 
tainty in estimating its age. Mannings & Sargent (1997) 
estimate a total stellar luminosity and effective temper- 
ature, = 53.6 Lq and T e fj = 10210-fT, which corre- 
sponds to a stellar age and mass of ^3-5 Myr and 2.5 M Q , 
using pre-main sequence tracks computed by D Antona & 
Mazzitelli (1994). Van den Ancker, de Winter, & Tjin A 
Djie (1998) have estimated a total stellar luminosity and 
effective temperature, — 47 L Q and T e // = 9500ii, 
which corresponds to a stellar age and mass of ^2 Myr 
and 2.4 Mq, using pre-main sequence tracks computed by 
Palla & Stahler (1993). AB Aur is an excellent source for 
high resolution mid-infrared imaging because it is bright 
and close by. It is one of four stars in the Hillenbrand 
et al. (1992) sample which possess a Q band magnitude 
brighter than -1.0 and lies at a distance less than 150 pc. 
All models for the dust distribution around AB Aur pre- 
sume that it is a single star. If AB Aur is a binary, then 



current models are incomplete. 

Millimeter aperture synthesis imaging of AB Aur has re- 
solved a rotating disk of 13 CO (J=l— >0), with a radius of 
^450 AU and at a position angle of 79°, around AB Aur, 
with velocities consistent with Keplerian rotation (Man- 
nings & Sargent 1997). Based upon the aspect ratio of 
the major:minor axes, Mannings & Sargent (1997) find 
an inclination of 76° for the disk. Previous mid-infrared 
imaging studies of AB Aur report spatially resolved ob- 
servations of the mid-infrared emission at 17.9 /im using 
the 5 m Hale telescope (Marsh et al. 1995). They find 
an elongated disk, with an East- West extent of 80±20 AU 
and a North-South extent of 25±20 AU and thus a posi- 
tion angle consistent with that observed for the CO disk. 
However, subsequent higher resolution near infrared inter- 
ferometry (Millan-Gabet et al. 1999; Millan-Gabet, Schlo- 
erb, & Traub 2001) and scattered light imaging of AB Aur 
(Grady et al. 1999) suggest a disk inclination of less than 
45°. 

The dust around AB Aur has been studied with a variety 
of techniques. Near infrared, interferometric observations 
at H and K' bands have resolved a ring of emission with 
a radius of 0.35 AU (Millan-Gabet et al. 1999). ISO 2- 
200 /im spectra have revealed the presence of PAH, FeO, 
and olivine emission features and Si-0 and O-Si-0 bending 
modes (van den Ancker et al. 2000). However, to date, 
no complete view of the dust around AB Aur has been 
constructed. Thus, we discuss models for the dust around 
AB Aur in the context of millimeter continuum fluxes, in- 
frared spectra, and near infrared visibilities. In addition, 
we compare models for the dust around AB Aur with our 
high resolution Keck 18.7 /im map. 

2. OBSERVATIONS 

Our data were obtained on 2000 February 20 (UT) and 
2000 August 9 (UT) at the Keck I telescope using the 
Long Wavelength Spectrometer (LWS) which was built by 
a team led by B. Jones and is described on the Keck web 
page. The LWS is a 128 x 128 SiAs BIB array with a pixel 
scale at the Keck telescope of 0.08" and a total field of view 
of 10.2"xl0.2". We used the "chop-nod" mode of observ- 
ing and two different filters: 11.2-12.2 /im and 18.2-19.2 
/xm. The seeing on 2000 August 9 was significantly better 
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than that on 2000 February 20; thus, we concentrate our 
analysis on the data obtained in August 2000. We used a 
Cet for flux and point spread function (PSF) calibrations 
and MWC 480 (an unresolved Herbig Ae star) and 7 And 
for additional PSF calibrations. The data were reduced at 
UCLA using standard LWS routines. 

We flux calibrate our February and August data using 
the result for a Tau and a Cet that F„(11.5/j,m) = 507.8 Jy 
and 174.5 Jy, respectively and F v (19.3fj,m) = 189.4 Jy and 
66.3 Jy (Gezari et al. 1987) extrapolated to our bands as- 
suming that F v oc v 2 between 10 and 20 /an. For AB Aur, 
we find F v (11.7fim) = 25±2 Jy and 19±2 Jy for February 
and August, respectively, and F u (18.7/j,m) = 31±2 Jy and 
17±2 Jy for February and August, respectively We con- 
servatively estimate the uncertainties associated with our 
measurements from the drift in the fluxes of all of our stan- 
dard stars throughout the night. Comparison of our fluxes 
and of published fluxes in the 10 /im and 20 /im bands 
(listed in Tables 2 and 3) suggest that AB Aur is variable 
at mid- infrared wavelengths. The 24±10% 11.7 /im vari- 
ability we observe is consistent with the 21.7±6.6% IRAS 
12 /zm 6 month variability reported by Prusti & Mitskevich 
(1994). The 46±8% 18.7 /im variability we observe is sig- 
nificantly larger than the 14.3±8.5% variability IRAS 25 
/im 6 month variability reported by Prusti & Mitskevich 
(1994). However, the absolute magnitudes of our 18.7 /im 
fluxes are similar to 17 /tni fluxes measured using ISO and 
the IRTF, 24.4 Jy (Thi ct al. 2001) and 22.7 Jy (Richter 
et al. 2002) respectively. 

The color temperature, inferred from our 11.7 /im and 
18.7 /im photometry, for our data agrees within the er- 
ror bars for our two epochs of observations. The flux ra- 
tio, F„(11.7 /im)/F„(18.7 urn), is 0.8±0.1 and 1.1±0.2 for 
Februrary 2000 and August 2000, respectively. Thus, we 
can not rule out the possibility that changes in intrinsic lu- 
minosity of the star generate the observed variability. AB 
Aur possesses strong variability in the Ha P Cygni profile 
and in measurements of the linear polarization of UBVRI. 
Based upon monitoring of these diagnostics in 1993 and 
1994, Beskrovnaya et al. (1995) suggest that circumstel- 
lar inhomogeneities may be responsible for the variability 
on hour to month long timescales. Prusti & Mitskevich 
(1994) suggest that dust formation in the stellar wind or 
dust on eccentric orbits may be reponsible for the mid- 
infrared variability. 

We generate images from our chop-nod sets by coadding 
chop pairs from both nod positions in each of our nod sets. 
The 18.7 fim PSF contains a ghost and diffuse asymmetric 
scattering while the 11.7 /im PSF is clean. Although the 
seeing was excellent, the PSFs varied somewhat during the 
night. We show line cuts through AB Aur and our 11.7 
/im standards in Figure la and through AB Aur and our 
18.7 /im standards perpendicular to and parallel to the 
scattered light in Figures lb and lc. The Airy rings can 
be seen in all of the 18.7 /im line cuts and half of the 11.7 
/im line cuts. Since there is some variation in the PSF 
at 18.7 /im, we compare our images of AB Aur to 3 PSF 
calibrators: a Cet, observed 57 minutes earlier, MWC 480 
(an unresolved Herbig Ae star), observed 11 minutes later, 
and 7 And, observed 2 hours, 26 minutes earlier than AB 
Aur. We use MWC 480 as our primary PSF because it 
was observed closest in time to AB Aur. We also use 7 



And because it was observed with the worst seeing during 
the evening. 

AB Aur is probably resolved at 18.7 /tm with an angular 
diameter of 1.2" at a surface brightness of 3.5 Jy arcsec~ 2 
(10% of the maximum surface brightness). An image of 
our calibrator a Cet at 18.7 /im is shown in Figure 2 with 
contours showing 1%, 2%, 3%, 5%, 10%, 25%, 50%, and 
99% of the maximum surface brightness. We display a Cet 
instead of MWC 480 because it is a significantly bright 
object, making the ghost and diffuse scattered light sig- 
nificantly easier to see. An image of AB Aur at 18.7 /im 
is shown in Figure 3 for comparison with the same con- 
tours as in the image of a Cet. The diameter of AB Aur 
at 10% of the maximum surface brightness is 50% larger 
than the same contour in the image of a Cet. In addition, 
the lowest surface brightness contours (1%, 2%, 3%, and 
5%) are more symmetric in the image of AB Aur than in 
the image of a Cet. Our image of AB Aur at 18.7 /tm 
does not exhibit the same asymmetry reported by Marsh 
et al. (1995) at 17.9 /tm. AB Aur appears unresolved at 
11.7 /zm, with a Full Width at Half Maximum (FWHM) 
of 0.30", consistent with a point source. 

3. THE CIRCUMSTELLAR DISK AND MORE? 

The total mass of dust around AB Aur has been esti- 
mated from millimeter contiuum measurements in an el- 
liptical beam with FWHM 3.5" x 5.6". At A = 2.7 mm, 
Mannings & Sargent (1997) find that AB Aur appears un- 
resolved and measured a flux, F„(2.7 mm) = 10.6 ± 0.4 
mJy which, in their models, corresponds to Mdust = 2.0 
x 10 29 g. 

We can estimate the minimum mass of warm dust con- 
tributing to the 18.7 /tm flux. If the dust is optically thin 
and the population has a single temperature, then 

F D 2 

M - — * m 

™ " B v {T gr ) Xv W 

(Jura et al. 1995) where \v 1S the dust absorption opacity 
and D* is the distance to the object. We estimate the tem- 
perature of the dust, T gr , from the ratio of our 11.7 /tm 
and 18.7 /im fluxes, assuming a single population of grains. 
If T gr = 250 K, F„(18.7 /im) = 17.2 Jy, and X v = 1000 
cm 2 /g (Ossenkopf, Henning, & Mathis 1992), then M warm 
= 1.2 x 10 25 g, substantially less than that inferred from 
millimeter continuum measurements. Therefore, a large 
amount of cold dust must exist in the system, presumably 
in a disk. 

4. COMPARISON WITH MODELS 

We use our data to evaluate the following four detailed 
models for the dust around AB Aur: (1) An optically thick, 
geometrically thin disk surrounded by an optically thin, 
spherically symmetric envelope (Vinkovic ct al. 2002). 
(2) Two spherically symmetric dust shells, whose grain 
properties (composition, size and spatial distribution) are 
inferred from ISO 2 - 200 /tm spectra (Bouwman et al. 

2000) . (3) A geometrically thin, flared disk in vertical hy- 
drostatic equilibrium (D'Alessio et al. 1998). (4) A flared 
disk with an optically thin surface layer of small grains and 
an infrared emitting rim (Dullemond, Domink, & Natta 

2001) . 
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We convolve these models for the predicted surface 
brightness of AB Aur at 18.7 /im with the PSF and com- 
pare them with our observations of AB Aur. Each model 
was interpolated to provide surface brightnesses on the 
same grid as the observations. Since line cuts through each 
of the models, convolved with a PSF, compared with the 
AB Aur data look similar independent of position angle, 
we only show line cuts along the North-South direction 
of each of the models, convolved with a PSF, compared 
with the AB Aur data. In Figure 4a, MWC 480, an un- 
resolved Herbig Ac star observed closest in time to AB 
Aur, was used as the PSF. In Figure 4b, 7 And, the most 
extended PSF observed during the night, was used. Line 
cuts through different position angles appear similar; thus, 
we show line cuts through only one orientation. In addi- 
tion, we calculate the reduced x 2 for the inner 1.2" x 
1.2" centered on AB Aur, using MWC 480 as the PSF. In 
both analyses, we scale the image of the model, convolved 
with the PSF, so that the maximum surface brightness is 
equal to the measured maximum surface brightness of AB 
Aur. Since AB Aur is located at a distance of 144 pc, 
the flux expected from the stellar photosphere at 18.7 pm, 
F v (18. 7 fim) — 0.032 Jy, is significantly less than the small- 
est measured excess flux, F v (18.7pm) = 17.2 Jy, and can 
be neglected. We calculate the reduced \ 2 f° r eac h model 
assuming that the uncertainty in the surface brightness is 
10%, approximately the uncertainty in the 18.7 pm pho- 
tometry. Finally, we also compare the models to measured 
near infrared interferometric visibilities, infrared spectra, 
and millimeter continuum measurements. 

4.1. A Spherically Symmetric Envelope? 

Vinkovic et al. (2002) have modeled the SED of AB Aur 
with a system which consists of a star with radius, i?* , and 
effective temperature, T#, surrounded by a geometrically 
thin and optically thick passive disk extending from i?„ to 
an outer radius Rdisk — 240 AU. The optically thick disk 
is assumed to have an inclination, i = 76°, consistent with 
the CO observations of Mannings & Sargent (1997). They 
also incorporate a spherical dust envelope with an inner 
radius, R su b = 99.3 i?* and an outer radius, 1000 R su b- 
At distances less than R su b, grains in the envelope are 
sublimated. The mass density and grain size distributions 
are given by 

2 



P(r) = Po 



R 



sub 



0-2Rsub 



n(a) = Aa 



-3.5 



(2) 
(3) 

where r is the distance between the star and the circum- 
stellar material and a is the grain radius. The dust grains 
are assumed to be an interstellar mixture of 53% silicates 
and 47% graphite by number (Draine & Lee 1984) with a 
minimum grain radius, a m in — 0.005 /im, and a maximum 
grain radius a max = 0.250 /im. They assume an optically 
thick, geometrically thin (flat) disk with an outer temper- 
ature of 25 K, a disk mass, Mdisk = 8.6 x 10 30 g, and an 
envelope mass, M env = 2.8 x 10 27 g, corresponding to an 
optical depth, ry = 0.44, consistent with the value mea- 
sured by van den Ancker et al. (2000) (Ay = 0.5). Their 
best fit model for the SED for AB Aur yields a total flux 
F v (18.7pm) = 25 Jy with 12 Jy from the disk and 13 Jy 
from the envelope. 



We compare the inner 1.2" x 1.2" of our 18.7 /im map 
of AB Aur with the model (Vinkovic et al. 2002). When 
the intensity map for this model is convolved with the 
PSF and scaled so that the maximum surface brightness 
is equal to that measured for AB Aur, the reduced x 2 = 
4.6. Their model agrees well with our data. 

However, the Vinkovic et al. (2002) model does have 
some limitations. They do not provide a physical expla- 
nation for the stability of the envelope. In addition, their 
model does not attempt to reproduce the detailed infrared 
spectra; a standard interstellar grain mixture is assumed. 
While the near infrared fluxes are fit, the measured visibil- 
ities are not considered. Finally, the model assumes a dust 
disk mass which is more than an order of magnitude higher 
than is inferred from millimeter continuum measurements. 

4.2. Modeling the Infrared Spectra 

Bouwman et al. (2000) infer dust composition, spatial 
distribution, and grain sizes by fitting ISO 2 - 200 pm 
spectra. For simplicity, they assume that the dust around 
AB Aur is optically thin and contained within spherically 
symmetric shells and neglect the presence of a circumstel- 
lar disk. They describe the density distribution, p(r), and 
the grain size distribution, n(a), with power laws. 
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where p a is the dust density at the inner edge of the dust 
shell, Ri n , and a min is the minimum dust grain size. They 
find that a combination of two populations of dust, a "hot" 
population at 1 - 11 AU, containing 5.2 x 10 24 g (m = 2.8) 
and a "cold" population at 28 - 175 AU, containing 6.6 x 
10 28 g (m = 2.0), reproduces the observed 2-200 ^m ISO 
spectrum well. Their model fit parameters, including grain 
composition and the mass fraction of dust in each species, 
are reproduced in Table 4. Bouwman et al. (2000) assume 
that the iron oxide grains are small (2ira/\ << 1) and do 
not infer a size distribution for this species. 

Since Bouwman et al. (2000) do not predict the surface 
brightness of their model, we estimate the surface bright- 
ness of optically thin, spherically symmetric dust shells as 
a function of angular offset on the sky, assuming the same 
"hot" and "cold" populations of dust with the same power 
law dependences for p(r) and n(a) and the same compo- 
sition. The specific intensity at an offset, b, from the star 
is calculated by integrating along the line of sight through 
the envelope, perpendicular to the sky. 



species 



/ Pi(r) / B„{Tg r s(x)}ni(a)Q abs sira 2 dada 

J X m i n J a,min,i 



^ _ _ (6) 

where r = \Jb 2 + x 2 is the distance between the star and 
the dust grains whose radiation is summed in the integral. 

For the cold grains, x max = y 175 AU 2 — b 2 and x m in = 

for b < 28 AU and x mm = \/28 AU 2 - b" 2 fo r b > 28 AU. 

For the warm grains,x maa; = y/Tl AU 2 - b 2 and 

for b < 1 AU and x min = y/l AU 2 - b 2 for b > 1 AU. 
The grain temperature for each species, T gr , is calculated 
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numerically, at a distance r from the star, assuming that 
the grains are in radiative equilibrium. 

( R \ 2 f°° f°° 

f-^J / Qabs(v)B v (T,)dv = 4/ Q abs (v)B u (T gr )dv 

(7) 

where T» and i?* are the stellar radius and temperature, 
D is the distance between the star and the dust grain. 
We calculate the absorption coeffecient for each species of 
dust, Qabs{v), and wavelengths less than 0.2 ^m follow- 
ing Bouwman et al. (2000). We calculate a total flux 
F„(18.7[j,m) = 45 Jy, with 20 Jy produced by the hot pop- 
ulation and 25 Jy produced by the cold population. 

The Bouwman et al. (2000) model does not fit our data 
particularly well. After scaling the surface brightness to 
that of AB Aur, convolving a 1.2" x 1.2" model of the 
dust emission with the PSF (MWC 480) yields a reduced 
X 2 = 8.6. 

Bouwman et al. (2000) use grain materials which repro- 
duce the ISO 2 - 200 /xm emission features well. They have 
some difficulty reproducing the observed near infrared visi- 
bilities. Their best fit model includes iron grains with radii 
0.01 - 0.1 /im to produce the near infrared emission close 
to the star. Iron is chosen as a material because it is one 
of the few materials which is likely to survive at tempera- 
tures ^1500 K and its emissivity is high at near infrared 
wavelengths and low at mid- and far infrared wavelengths; 
there is no direct evidence for the presence of metallic iron 
from emission line features. The Bouwman et al. (2000) 
model is inconsistent with the measured visibility curves. 
A single grain population of metallic iron grains with ra- 
dius, a = 0.32 fim, and minimum distance, R m i n — 0.35 
AU fits the visibilities the best and reproduces the mea- 
sured near infrared size of AB Aur. Bouwman et al. (2000) 
neglect the presence of a circumstellar disk altogether. 

4.3. A Flared Disk in Hydrostatic Equilibrium? 

Another possibility is that the flat SED is produced 
by a flared disk. Flared disks intercept more stellar flux 
than geometrically flat disks, allowing surface grains to be 
warmed to higher temperatures. Flared disks have been 
observed around the young stellar objects HH 30 and HK 
Tau/c in scattered light using HST (Burrows et al. 1996; 
Stapclfeldt 1998). Simple semi- analytical models of flared 
disks with a layer of small (2ira/\ « 1), optically thin, 
surface grains can reproduce the SEDs of some T Tauri 
and Herbig Ae stars (Chiang & Goldreich 1997; Chiang et 
al. 2001). Detailed numerical modeling has shown that 
stellar radiation is capable of heating the outer regions of 
the disk so that they become optically thin and vertically 
isothermal (D'Alessio et al. 1998). D'Alcssio, Calvet, & 
Hartmann (2001) successfully model the SEDs and images 
of edge-on disks like HH30 and HK Tau/c. 

D'Alessio et al. (2002) model the intensity of the dust 
around AB Aur at 18.7 /im assuming that the central star 
has mass, = 2.4 M Q , radius, i?* = 2.5 Rq, and ef- 
fective temperature = 9500 K and is surrounded by 
a disk with inner radius, R m in = 3 i?* and an outer ra- 
dius, Rdisk = 200 AU. They assume a disk inclination, i 
= 75.5°, and a disk dust mass, Mdust = 5.6 x 10 29 g, with 
a gas:dust ratio of 100:1. They assume that the disk is 
composed of silicates, trolite, ice, and organics with mass 
fractional abundances (relative to the gas mass) of 0.0034, 



0.000768, 0.0056, and 0.0041 respectively (D'Alessio et al. 
2001). The grains are assumed to have a size distribution 
described by a power law, 

n{a) = Aa- 3 - 5 (8) 

with a minimum grain radius, a m ; m = 0.005 /im and a 
maximum grain radius, a max = 1 mm. The assumption 
of grain growth in this model changes the opacity of the 
dust grains at small and large wavelengths and therefore 
the temperature distribution of the disk. Larger grain sizes 
decrease the fraction of stellar radiation intercepted by the 
disk resulting in a colder disk. The model assumes an ac- 
cretion rate, M = 1.0 x 10~ 8 M Q /year, corresponding to 
a accretion luminosity of 2 Lq . 

The flared disk model does not fit our data particularly 
well. It predicts too little 18.7 fim flux, 5 Jy, compared 
with all of the 20 /im band flux measurements. This may 
partly be the result of the high assumed inclination, i = 
75.5°. Recent near infrared interferometry (Millan-Gabet 
et al. 1999) and scattered light imaging (Grady et al. 
1999) suggests that the inclination of the AB Aur disk is 
probably less than 45°. Since the luminosity of the disk 
depends directly on the cosine of the inclination, the esti- 
mated flux of the D'Alessio disk might be as much as 2.8 
times higher if it is calculated for an inclination of 45°. 
When the maximum surface brightness is scaled to that 
of AB Aur, convolving a model of the inner 1.2" x 1.2", 
with the most likely PSF (MWC 480) yields a reduced \ 2 
= 9.2. However, convolving the model with the worst PSF 
(7 And), yields a model which appears consistent with the 
observations. 

The D'Alessio et al. (2002) model has some limitations. 
They assume a Pollack et al. (1994) grain mixture and ig- 
nore detailed grain composition information from infrared 
spectra. They do not fit the observed near infrared vis- 
ibilities. However, their assumed circumstellar dust disk 
mass is an approximate agreement with that inferred from 
millimeter continuum measurements. 

4.4. A Flared Disk with a Hot Inner Rim? 

The simpled flared disk model with a surface layer of 
small, optically thin grains predicts too little disk emis- 
sion for AB Aur at near infrared wavelengths. Dullcmond 
et al. (2001) have refined the Chiang k, Goldreich (1997) 
flared disk model to allow for direct stellar heating of dust 
grains at the inner disk radius. In their model, grains sub- 
lime when their temperatures reach T r i m (= T sub = 1500 
K). The direct stellar radiation on the inner rim surface 
warms the grains effectively, creating the observed near 
infrared excess, puffing up the inner disk, and casting a 
shadow on adjacent portions of the disk. Dullcmond et 
al. (2001) model the AB Aur system with a star, with M* 
= 2.4 M , i?» = 2.5 Rq, T eff = 9520 K, (L, = 47 Lq) 
surrounded by a disk with an inner rim R r i m = 0.52 AU 
away from the star and an outer radius Rdisk = 400 AU. 
The vertical height of the inner rim is H r i m = 0.19 R r i m - 
In modeling the SED, the circumstellar disk is assumed to 
have an inclination, i = 65°, less than that inferred from 
the CO observations of Mannings & Sargent (1997), more 
consistent with the 45° inclination inferred from scattered 
light observations of Grady et al. (1999). However, the 
model of intensity as a function of position was calculated 
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assuming that the AB Aur disk was face-on. Dullemond 
et al. (2001) assume a disk surface density 

E = E »(rSu)~ 2 W 

where S D = 10 4 g/cm 2 with a total mass of 4.0 x 10 29 g. 
The gas:dust ratio is assumed be 100:1. The dust grains 
are assumed to be a mixture of 95% astronomical silcate 
(Draine & Lee 1984) and 5% amorphous carbon (Ivezic 
et al. 1997). The rim and disk produce an 18.7 ^m flux, 
^(18.7/xm) = 49 Jy. 

The Dullemond et al. (2001) model does not agree with 
our mid- infrared observations. The reduced \ 2 f° r the 
model convolved with the PSF (MWC 480), in the inner 
1.2" x 1.2", is 84, when the convolved image is scaled to 
have the same maximum surface brightness as AB Aur. 
Since the rim produces an order of magnitude less emis- 
sion than the disk at 18.7 /im, the Dullemond et al. (2001) 
model is essentially a Chiang & Goldreich (1997) flared 
disk with an inner radius of 8 AU. If the disk were to 
extend to an inner radius of 0.5 AU, 8% more 18.7 ^im 
emission would be produced (Dullemond et al. 2001) in 
the central pixel, only slightly affecting our normalization. 
Thus, the Chiang & Goldreich (1997) model similarly does 
not agree with our mid-infrared observations. 

Dullemond et al. (2001) do fit other observational data 
well. They reproduce the measured near infrared fluxes 
and visibilities by placing the hot inner rim at the grain 
sublimation distance. Their assumed disk mass agrees 
with the dust mass, inferred from millimeter continuum 
measurements. They use a silicate/carbon grain mixture 
and do not match the infrared spectrum. 

5. DISCUSSION 

Models for the dust around AB Aur must fit a suite of 
multi wavelength observations if they are to be successful. 
In Table 5, we summarize the current ability of models to 
fit various observational data sets. In Table 6, we com- 
pare the predicted 18.7 fim fluxes, dust disk masses, and 
reduced % 2 s of each model to those inferred from obser- 
vations. By examining the line cuts shown in Figures 4a 
and 4b, we see that the Dullemond et al. (2001) model 



for the dust intensity predicts too much flux at distances 
>0.24" (35 AU) from the star. If MWC 480 more accu- 
rately represents the PSF, then the Bouwman et al. (2000) 
model predicts too much emission and the D'Alessio (2002) 
model predicts too little emission at distances >0.40" (58 
AU) from the star. In this case, the Vinkovic et al. (2002) 
model is favoured. Flared disk models are not able to re- 
produce both the compact size of AB Aur and the strong 
mid-infrared flux. However, there is some uncertainty in 
the seeing at the time that AB Aur was observed; thus, it 
is not possible to rule out the flared disk model. Limita- 
tions exist with all the models, as described above. Future 
models of AB Aur must consider the grain composition, if 
they are to accurately model the SED and dust emission. 

6. CONCLUSIONS 

We have obtained high resolution mid-infrared images 
of AB Aur at 11.7 ^m and 18.7 /im using the LWS on the 
Keck I telescope. 

1. AB Aur is probably resolved at 18.7 /tm, with an 
angular diameter of 1.2" at a surface brightness of 3.5 Jy 
arcsec~ 2 , and unresolved at 11.7 /um, with a FWHM = 
0.30". 

2. Only ^10 25 g of dust is required to produce the 18.7 
/im emission. Millimeter interferometry indicates the pres- 
ence of 2 x 10 29 g of dust. Thus, a large amount of cold 
dust must exist in the system, presumably in a disk. 

3. The spatial distribution of the 18.7 ^m emission is 
most consistent with radiation from an optically thick, ge- 
ometrically thin disk embedded in a spherically symmetric 
dust envelope (Vinkovic et al. 2002); however, the pres- 
ence of a flared disk can not be excluded. None of the 
models yet reproduces all of the available observations. 
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Table 1 
AB Aur Properties 

Quantity Adopted Value Reference 



Spectral Type 


AOVc+sh 


1 


Distance 


144 pc 


2 


Effective Temperature (T e ff) 


9500-10210 K 


2, 5 


Stellar Radius (R„) 


2.3-2.5 R© 


2, 5 


Stellar Luminosity (L*) 


47-53.6 L Q 


2, 5 


Stellar Mass (M*) 


2.4-2.5 M Q 


2, 5 


Rotational Velocity (usini) 


80 km/sec 


3 


Fractional Dust Luminosity 


0.44 


4 


(W£.) 






Estimated Age 


2-5 Myr 


2, 5 



References. — (1) Bohm & Catala (1993); (2) van den An- 
cker et al. (2000); (3) Bohm & Catala (1995); (4) Natta et al. 
(2001); (5) Mannings & Sargent (1997) 
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Table 2 



Measured 10 Fluxes 



Date 


Wavelength 


Flux 


Reference 


Observed 


(A*m) 


(Jy) 






10.1 


25.2 


Hillenbrand et al. (1992) 


1993 Oct 23 


11.7 


20.6±0.6 


Marsh et al. (1995) 


1994 Nov 14 


11.7 


23.0±0.7 


Marsh et al. (1995) 


2000 Feb 20 


11.7 


25±2 


this paper 


2000 Aug 9 


11.7 


19±2 


this paper 




12 


30.4±1.8 


Prusti & Mitskevich (1994) 




12 


23.8±1.4 


Prusti & Mitskevich (1994) 




12 


18.7±1.3 


Richter et al. (2002) 



Table 3 
Measured 20 fiM Fluxes 



Date 


Wavelength 


Flux 


Reference 


Observed 


(/mi) 


(Jy) 






17 


24.4 


Thi et al. (2001) 




17 


22.7±2.3 


Richter et al. (2002) 


1994 Nov 14 


17.9 


34.7±5.2 


Marsh et al. (1995) 


2000 Feb 20 


18.7 


31±3 


this paper 


2000 Aug 9 


18.7 


17±2 


this paper 




20.25 


27.2 


Hillenbrand et al. (1992) 




25 


49.6 


Prusti & Mitskevich (1994) 




25 


42.5 


Prusti & Mitskevich (1994) 



Table 4 



Grain Properties 





\A/ enrol on erf" n 
VV ctvclcllg, 111 




±\.L,1L,I L.I1L-L, 




Observed 


( um) 


(Jv) 






Vinkovic et al. (2002): 












Mf rac 


a (/jm) 






Astronomical Silicate 


0.63 


0.005-0.250 






Carbon 


0.37 


0.005-0.250 






Bouwman et al. (2000): 












Hot M frac 


Hot a (/Um) 


Cold M frac 


Cold a (fim) 


Olivine 


0.71 


0.01-5.0 


0.74 


0.01-126 


Carbon 


0.16 


0.01-2.0 


0.11 


0.01-32 


Water Ice 






0.15 


0.1-40.0 


Iron 


0.06 


0.01-1.0 






Iron Oxide 


0.07 




9xl0~ 4 




D'Alcssio (2002): 












Mf rac 


a (/an) 






Silicate 


0.25 


0.005-1000 






Trolite 


0.05 


0.005-1000 






Ice 


0.40 


0.005-1000 






Organics 


0.30 


0.005-1000 






Dullemond et al. (2001): 












Mfrac 


a (/mi) 






Astronomical Silicate 


0.95 


0.005-0.250 






Amorphous Carbon 


0.05 


0.16 







Table 5 

Comparing the Models with Observations 





Infrared 


18.7 /j,m 


Near infrared 




Spectrum 


Image 


Visibilities 


Vinkovic et al. (2002) 


Omitted 


Yes 


Omitted 


Bouwman et al. (2000) 


Yes 


No 


No 


D'Alessio (2002) 


Omitted 


Possibly 


Omitted 


Dullemond et al. (2001) 


Omitted 


No 


Yes 



Table 6 
Model Properties 





F u (18.7 /urn) 


Dust Mass 


Reduced 




(Jy) 


(g) 


x 2 


Observed Quantities 


31, 17 


2 x 10 29 t 




Predicted Quantities: 








Vinkovic et al. (2002) 


25 


8.6 x 10 30 


4.6 


Bouwman et al. (2000) 


45 


6.6 x 10 28 t 


8.6 


D'Alessio (2002) 


5 


5.6 x 10 29 


9.2 


Dullemond et al. (2001) 


49 


4.2 x 10 29 


84 



^Disk mass estimated from 2.7 mm interferometry (Mannings & Sar- 
gent 1997) 

-^Dust in the Bouwman et al. (2000) model is located in spherically 
symmetric dust shells and not in a circumstellar disk. 
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Offset (arcsec) Offset (arcsec) 

Fig. 1. — (a) Line cuts through AB Aur (black) and our 11.7 fim standard stars 7 Dra (red) and a Cet (green), (b) Line cuts through AB 
Aur and our 18.7 (im standard stars 7 Dra (red), a Cet (green), and MWC 480 (blue) made perpendicular to the direction of the asymmetric 
scattered light, (c) Line cuts through AB Aur and our 18.7 fim standard stars 7 Dra (red), a Cet (green), and MWC 480 (blue) made parallel 
to the direction of the asymmetric scattered light. 
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-4 -2 2 4 

OFFSET (ARCSEC) 



Fig. 2. — The 18.7 /an image of a Cet, our standard star, rotated so that the ghost appears in the same position as in the image of AB 
Aur. The contours show 1%, 2%, 3%, 5%, 10%, 25%, 50%, 75%, and 99% of the maximum surface brightness. 
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Fig. 3. — The 18.7 /im image of AB Aur. North is up and East is to the left. The color bar shows the intensity key (Jy arcsec 2 ). The 
contours show 1%, 2%, 3%, 5%, 10%, 25%, 50%, 75%, and 99% of the maximum surface brightness. 



-2-10 1 2-2-10 1 2 

Offset (arcsec) Offset (arcsec) 

Fig. 4. — (a) Line cuts North-South through AB Aur at 18.7 /im and each of the models convolved with MWC 480, the PSF observed closed 
in time to AB Aur. (b) Line cuts North-South through AB Aur at 18.7 (im and each of the models convolved with 7 Dra, the standard star 
observed with the worst seeing. The blue, red, green, and cyan lines (from bottom to top )show the D'Alessio (2002), Vinkovic et al. (2002), 
Bouwman et al. (2000), &: Dullcmond et al. (2001) models respectively. 



